Strongly enhanced thermal transport in a lightly doped Mott insulator at low 
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We show how a lightly doped Mott insulator has hugely enhanced electronic thermal transport at 
low temperature. It displays universal behavior independent of the interaction strength when the 
carriers can be treated as nondegenerate fermions and a nonuniversal "crossover" region where the 
Lorenz number grows to large values, while still maintaining a large thermoelectric figure-of- merit. 
The electron dynamics are described by the Falicov-Kimball model which is solved for arbitrary 
large on-site correlation with a dynamical mean-field theory algorithm on a Bethe lattice. We show 
how these results are generic for lightly doped Mott insulators as long as the renormalized Fermi 
liquid scale is pushed to very low temperature and the system is not magnetically ordered. 
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Introduction Thermoelectric materials are attracting 
significant attention, because of their potential for vari- 
ous power generation or refrigeration applications which 
involve so-called green technologies. The mass applica- 
tion of the thermoelectric devices is hampered by their 
low efficiency and the aim of current research is to pro- 
duce materials with better thermoelectric conversion ef- 
ficiency. At high temperatures, the efforts are directed 
towards nanostructured semiconductors [l[ with reduced 
thermal conductivity. At low temperatures, the focus is 
on the materials with strongly correlated electrons, like 
Kondo insulators and systems with a Mott-Hubbard gap 
to enhance the thermopower in metallic systems. 

Here, we present a theory for the charge and thermal 
transport in a slightly doped Mott insulator that is de- 
scribed by the Falicov-Kimball model. The exact solution 
is obtained from dynamical mean-field theory (DMFT) 
and it shows that for large correlation and small doping 
the figure-of-merit is unusually large. This behavior is 
expected to also hold in more general doped Mott insu- 
lators, like those described by the Hubbard model if they 
are sufficiently frustrated away from magnetic order and 
if the operating temperature lies well above the renormal- 
ized Fermi temperature. Even the presence of magnetic 
order is unlikely to strongly modify this effect, since the 
magnetic order rarely has a large effect on charge and 
heat transport in strongly correlated materials [2|. The 
numerics, however, requires quite high precision to de- 
termine the transport for such low doping values, which 
is why this problem can only be solved within the Falicov- 
Kimball model with current state-of-the-art methods. 

Theoretical description To discuss this thermoelectric 
phenomena, we use transport equations which express 
the charge and the internal energy current densities, J(x) 
and Jg(x), in terms of the generalized forces Q. The co- 
efficients of these generalized forces, Nij(T), are given by 
various current-current correlation functions which have 
to be calculated for the model at hand. The electrical 



conductivity er, the Seebeck coefficient (or thermopower) 
a, and the thermal conductivity then obtained as 
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where D = N22/N11 — N^/N^ gives the effective Lorenz 
number £(T) = K e joT = {e/k B ) 2 D{T)/T 2 . The dimen- 
sionless figure-of-merit of a particular thermoelectric ma- 
terial is ZT — a 2 o~T / n, where k = n e + n p h is the overall 
thermal conductivity due to the electronic and the lattice 
degrees of freedom. The electronic figure-of-merit can be 
expressed as ZT = a 2 / C and an efficient thermoelectric 
material has ZT > 1. We do not consider the phonon 
contribution to the thermal conductivity further here. 

We next show how to exactly evaluate the transport co- 
efficients for the spinless Falicov-Kimball model [i ] with a 
large on-site Coulomb interaction. The Hamiltonian on 
an infinite-coordination (Z — > 00 with Z the coordina- 
tion number in this equation only and not to be confused 
with the figure of merit ZT) Bethe lattice is 

n = -^Y,( c h+ h - c -y^+^)Y, c i c i +u Y, c i c i Wi ■ 
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where cj (cj creates (destroys) a conduction electron at 
site i, t* is the renormalized hopping, [i is the chemi- 
cal potential shifted so that \i = corresponds to half- 
filling, U is the interaction between localized and itiner- 
ant electrons, and Wi is a classical variable equal to if 
a localized electron is not at site i and equal to 1 if the 
localized electron is at site i. The average filling of the 
localized electrons is fixed at 1/2, so that the system is 
in a Mott-insulating state when U is larger than 2t* and 
fj, = 0. The sum over i and j in the kinetic-energy term 
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FIG. 1: (Color online) The density of states plotted as a 
function of energy for U = 2.2 (full line), U — 4 (dashed 
line), and (7 = 5 (dashed-dotted line). The lower and the 
upper gap-edges are at u)-(U) and cj+(U), respectively. 



is over nearest neighbor pairs and h.c. denotes the her- 
mitian conjugate. This model is appropriate for doped 
Mott insulators where the dopants provide carriers with- 
out producing impurity levels in the gap, as is often seen 
in many semiconductors. 

To calculate the thermal transport, we need to evaluate 
a series of Kubo formulas for the relevant transport co- 
efficients. This is greatly simplified because the Falicov- 
Kimball model satisfies the Jonson-Mahan theorem 
Thus, the transport coefficients can be written as 
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with the transport function given by 



A tr (u,) = —t I de Po (e)(4^ - ^)[Im G{e,u)Y , (6) 



where the noninteracting density of states is 
1 
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and G(e,Lu) is the "band-energy-dependent" Green's 
function (e is the noninteracting band energy) which we 
calculate within DMFT 0. (One should note that there 
is no momentum on a Bethe lattice, and instead, one 
should think of the band-energy as analogous to the inde- 
pendent variable of momentum.) The integral in Eq. 
can be performed exactly, resulting in 



Me) = ^ [ G („)] LLiL_ . 
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The local Green's function, G(lj) = J dep (e)G(e,uj) sat- 
isfies a simple cubic equation [9|, which allows for the 



FIG. 2: (Color online) Chemical potential plotted as a func- 
tion of temperature. The black, red, purple, green, and blue 
curves are calculated for U =2.2, 2.5, 3, 4, and 5, respectively. 
The full, dashed, and dashed-dotted curves correspond to the 
concentrations of n c = 10 -6 , 10 -5 , and 10 -4 electrons above 
half filling. The characteristic temperature T M , obtained for 
U — 4, is indicated by the full, dashed and dashed-dotted ar- 
rows. The inset show the low-temperature behavior obtained 
for U =4. 



numerics to be carried out to high precision (i.e., there is 
no self-consistent iterative algorithm needed to solve the 
problem). In this equation, and in the following, we set 
t* = 1 as our energy unit. Note that because the density 
of states at the band edge grows like the square root of 
frequency, the transport density of states is linear near 
the band edge as one goes into the band and vanishes as 
one goes into the gap. Similar results are obtained for 
a three-dimensional cubic lattice (within DMFT) , as has 
been checked for a few points in parameter space. 

Numerical results The energy dependence of the 
renormalized density of states [p(cj) = — ImG(w)/7r] of a 
Mott insulator is shown in Fig. [T] for several values of U, 
with /i = 0. The density of states for the Falicov-Kimball 
model does not depend on temperature, and since the cu- 
bic equation for G(lo) depends only on uj + fi, U, and the 
average filling of the localized electrons, the density of 
states for a lightly doped Mott insulator is identical to 
that of the Mott insulator, except the origin is shifted, 
as the chemical potential changes. For U > 2, the DOS 
is split into a lower and upper Hubbard band of width 
W and with a separation between the maxima approxi- 
mately equal to U. The shape of the Hubbard bands is 
nearly independent of U (for U > 2) and only the gap, 
which extends from w_(C7) to lj+(U), increases with U. 

The temperature dependence of the chemical poten- 
tial, obtained for various values of U and several (small) 
values of the upper Hubbard band filling n c , is shown 
in Fig. (The total conduction electron filling is just 
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1/2 + n c .) At zero doping, the system is electron- hole 
symmetric and the chemical potential is in the middle 
of the gap, fi(T) = 0. At finite electron doping, the 
zero-temperature chemical potential is in the upper Hub- 
bard band, just above the band edge, as can be seen 
by comparing /j,(0) in Fig. [5] with oj + (U) in Fig. [TJ For 
n c < 1, the values of /i(0) are approximately given by 
u+(U) rj (U — W)/2. The low-temperature behavior 
of fi(T) is demonstrated in the inset of Fig. [2] The 
data show that after an initial parabolic decrease, /x(T) 
is nearly linear up to a characteristic temperature T M . 
For T > Tp, the decrease of fiiT) slows down, because 
excitations across the gap restore the electron-hole sym- 
metry and the chemical potential approaches the high- 
temperature limit, jLt(T) — > 0. For a given U, the charac- 
teristic temperature T M increases with n c . This is shown 
in Fig. [2] by the full, dashed, and dashed-dotted arrows 
which indicate 7), obtained for U = 4 and n c = 10~ 6 , 
10 -5 , and 10 -4 , respectively. 

The observed behavior follows from the fact that n(T) 
is close to the bottom of the upper Hubbard band at 
low temperatures and that the excitations across the gap 
establish a symmetric state at high temperatures, where 
fi s» 0. Thus, for a given n c , an increase of U shifts /x(0) 
and T M to higher values and translates n(T) upwards, as 
shown in Fig. [21 If we increase n c keeping U constant, a 
higher temperature is needed to restore the electron-hole 
symmetry, so that T M increases with n c . 

To calculate the transport integrals in Eq. ([5]) we adjust 
/x(T) to yield the target filling of particles and measure 
the energy at each temperature with respect to /u(T). 
Since neither p(to) nor A tr (ui) change their shape, the 
transport coefficients in Eq. ([5]) are easy to compute. At 
high temperatures, T > 2),(J7), the behavior is universal, 
i.e., the transport coefficients are independent of n c , for 
a given U. This is revealed most clearly by the effective 
Lorenz number, C(T), plotted in the uppermost panel 
of Fig. [3] The curves C(T) obtained for a given U and 
various n c merge for T > T^, due to the fact that at 
such high temperatures [i ~ and all the systems we are 
concerned with become (nearly) electron-hole symmetric. 

At the lowest temperatures, T < T\, the Lorenz 
number is given by the usual value, C(T) ~ Cq = 
(7r 2 /3)(/cB/e) 2 , expected from the Wiedemann-Franz law 
for degenerate fermions. The characteristic temperature 
T\ depends on the distance of the T = chemical poten- 
tial from the upper Hubbard band edge //(0) — u> + (U). 
For very low doping, this temperature is physically irrel- 
evant and we do not show it on the figure. Another type 
of universality sets in when the chemical potential is in 
the gap and the lower Hubbard band starts affecting the 
transport. Figure [3] shows that, for a given U, the Lorenz 
number becomes independent of n c for a wide range of 
temperature T\ -C T < Tq, where Tq depends on the size 
of the Mott-Hubbard gap. Here, the Fermi gas is nonde- 
generate and the Wiedemann-Franz law holds with the 
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FIG. 3: (Color online) Lorenz number in units of [ks/e] 2 , 
thermopower in units of [fes/e], and thermoelectric figure of 
merit ZT are plotted as a function of temperature for the 
same parameters as in Fig. [2] 



smaller classical value Co — 2(fcB/e) 2 . 

Curiously, for T\ -C T < Tq, the thermopower and 
the figure of merit also assume universal forms, indepen- 
dent of U, which can be seen by the overlapping curves 
in Fig. [3l This low-temperature universality arises be- 
cause the shape of A tr (w) around u> = \i does not depend 
strongly on U . The universality is lost for T > Tq, when 
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the excitations across the gap become important. 

For To < T < T M , there is a crossover from the low- to 
the high-temperature regime. Here, the transport coeffi- 
cients vary in a non-universal way, which is also revealed 
most clearly by the effective Lorenz number. It grows 
rapidly for T > To(n c ), attains a maximum at about 
T ~ T M , and then decays slowly. Well above T M , where 
n and a assume their high-temperature universal forms, 
C(T) is independent of n c . 

We now consider in more detail the results for ther- 
mopower and the figure-of-merit, shown in Fig. [3J At 
very low temperatures, the slope of a(T) is determined 
by the logarithmic derivative din A tr /duj\ u= ^ which is 
large for low doping, since /i(T) is close to the band 
edge. Here, a(T) grows rapidly and a/T increases as 
n c decreases but, for a given n c , it does not depend on 
U. At higher temperatures, the presence of the Mott- 
Hubbard gap leads to novel features. Unlike in ordinary 
semiconductors, a(T) continues to grow when the tem- 
perature renormalization brings fJ>(T) below the bottom 
of the upper Hubbard band. When that happens, the 
states below the chemical potential do not contribute to 
N12 and a(T) grows to very large values for large U and 
T < T . Eventually, for T < T ~ T M , the system be- 
comes nearly electron- hole symmetric and a(T) drops to 
very small values. 

The figure-of-merit becomes very large in the proxim- 
ity of the metal-insulator transition. For U > 20 and 
doping which gives less than 10~ 19 conduction electrons 
(n c w 1(T 4 ), we find ZT > 200! The maximum of ZT 
occurs at a temperature at which the Wiedeman-Franz 
law holds for nondegenerate fermions and the thermal 
conductivity is very low, so that the thermal current due 
to any additional degrees of freedom will greatly reduce 
the thermoelectric efficiency. However, the temperature 
range in which ZT > 1 increases with U and n c . For 
larger doping, a relatively large ZT is found in a broad 
temperature range, even for moderate U. By tuning the 
concentration of conduction electrons, one can find a tem- 
perature window in which Cq <C C(T) and ZT ^> 1. 
Here, the electronic thermal conductivity is not small 
and the large ZT values should be physically relevant, if 
one were to include phonon effects to the thermal con- 
ductivity. This occurs, typically for temperatures of the 
order of T M . 

Summary and conclusions We presented a theory for 
the charge and heat transport of a "bad metal" described 
by the Falicov-Kimball model. Using DMFT, we show 
that the transport coefficients of a lightly doped Mott 
insulator exhibit universal features at low temperatures, 
Ti < T < T , and at high temperatures, T > T M . In the 
crossover region, To < T < T^, the effective Lorenz num- 
ber becomes very large and the Wiedemann- Franz law is 
not obeyed. The thermopower and the figure of merit 
grow to very large values close to T To. In this region, 
the system might be optimal for use in thermoelectric 



devices. 

The low-temperature thermoelectric properties are en- 
hanced because the renormalized transport density of 
states of a lightly doped Mott insulator is very asymmet- 
ric close to the gap. As the temperature increases, the 
chemical potential drops below the bottom of the con- 
duction band, such that the asymmetry of the electron 
and hole currents becomes even more pronounced, and 
the thermopower is further enhanced. At about T ~ T M , 
the chemical potential approaches the center of the gap, 
the system nearly acquires electron-hole symmetry, and 
the thermopower drops to small values. We could not 
obtain similar effects with non-interacting electrons or in 
a Fermi liquid. 

The ZT of a strongly correlated system with just a 
few electrons in the upper Hubbard band turns out to 
be surprisingly large. This finding is not restricted to a 
Bethe lattice; large ZT is obtained for a three-d cubic 
lattice as well. We point out that the maximum of ZT 
occurs at temperatures at which the Wiedemann-Franz 
law holds for nondegenerate fermions and the thermal 
conductivity is low, so that the thermoelectric efficiency 
of a real system described by our model would be greatly 
affected by other degrees of freedom that can transport 
heat. However, for large U and moderate n c , there is 
a broad temperature interval, T) < T < T M , in which 
ZT is moderately large, even though the electronic ther- 
mal conductivity and C are not small and in this regime, 
phonons will not affect the ZT values as much. 

One might ask whether such features are generic for 
all Mott insulators, or specific to the Falicov-Kimball 
model. It turns out, the Hubbard model will share this 
behavior for small enough doping and large U. This is 
because in this regime the renormalized Fermi tempera- 
ture is strongly reduced towards zero, so a Fermi liquid 
won't form in this interesting regime where the thermo- 
electic effects can be so large, and the Hubbard model 
acts similar to the Falicov-Kimbal model. (It is still pos- 
sible that a too large phonon thermal conductivity, that 
magnetic order which modifies the transport density of 
states, or localization effects due to disorder in lightly 
doped systems will reduce this effect, but there should 
be a regime where it will still be able to be seen ficj.) 
We hope that this work will inspire experimental groups 
to look for this kind of phenomena as a new route towards 
low-temperature thermoelectric refrigeration. 
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